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Abstract 

We present a derivation of the bosonic contribution to the thermodynamical 
potential of four fermion models by means of a l/iV c -expansion of the functional 
integral defining the partition function. This expansion turns out to be particularly 
useful to correct the mean field approximation expecially at low temperatures, where 
the relevant degrees of freedom are low-mass bosonic excitations (pseudogoldstones). 



1 Introduction 

Four fermion models, like the NJL M and the Gross- Neveu model P| , have been widely 
studied, expecially in connection with chiral symmetry breaking and restoration at finite 
densities and temperatures. The thermodynamics of these models may offer useful hints 
in understanding the properties of thermal QCD, expecially those related to the chiral 
symmetry. Such a knowledge is expected to be important to interpret heavy-ion collisions 
at relativistic energies and for the description of early phases of cosmology. 

Usually such studies are carried out in the mean field approximation, for instance by 
using the Dolan-Jackiw formalism || of imaginary-time Green functions. 

* Partially supported by the Swiss National Foundation 
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A review of the applications of the NJL model to hadronic and nuclear physics can be 
found in the report by Vogl and Wise 0, to which we also refer for the vast literature 
on the subject. A calculation of the thermodynamical potential for the NJL model in the 
grand canonical ensemble, in the mean field approximation, can be found for instance in 
ref-i- 

The mean field approximation neglects fluctuations around it, and it excludes in such 
a way possible bosonic contributions to the thermodynamical potential. On the other 
hand such fluctuations, as they stand out from possible resonant or non resonant bosonic 
scattering states, can be expected to be of great importance in the physical problem. For 
instance the pions, which are pseudogoldstones related to spontaneous chiral symmetry 
breaking, are expected to dominate the low energy properties of the system, as they corre- 
spond to the lowest mass bosonic excitations. We recall that in non-relativistic statistical 
mechanics a resembling situation is the calculation of the second virial coefficient, which 
is related to time-delay in the collision, and, in particular, bound states (two molecules 
combining to form a new molecule) lead to a modification of the system pressure ||. 

Unfortunately, even in a simple field theoretic model, such as NJL, exact analytical 
calculations appear unfeasible and one has to resort to approximate methods, such as 
large-N expansion. Such type of expansion has been used since a long time in statistical 
field theory. Well-known applications are to the O(N) symmetric field theory [|7| and to 
the Gross- Neveu model ||. 

A calculation of fluctuations has been recently carried out || for the NJL model by 
evaluating the thermodynamical potential through the method of coupling constant inte- 
gration, as developed in the book of Fetter and Walecka jlCl , by including terms to order 
1/N C , where N c is the number of colors (N c = 3 in QCD). The 1/N C contributions are 
taken into account through diagrammatic method. The use of 1/N C expansion goes back 



to Witten 11] and a review is given in the Erice lectures by Coleman |12 



The use of diagrammatic methods, besides its cumbersome features, has also a delicate 
aspect, as it may easily destroy general symmetries of the system by improper selection of 
the set of diagrams |13|- It is therefore strongly preferable to recur to procedures which 
automatically guarantee conservation of the system symmetries. 

In this note we shall develop a method based on the integral functional formulation of 
field thermodynamics which leads to straightforward calculation of the 1/N C expansion of 
the partition function. Within this method conservation of symmetries is consistently and 
naturally guaranteed. 
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2 1/N C expansion 

Let us consider a model in D = d + 1 dimensions defined by the following four fermion 
lagrangian 0. ^ 



C = if, (id - M) if) + 9 — [(^f - {hsTif)) 2 



(1) 



where M is the bare fermion mass, Ti are the matrices of the fundamental representation 
of the flavor group SU(Nf), and ip consists of iV c separate iVy-dimensional representations 
of SU(Nf). 

The generating functional (evaluated at vanishing sources) is 



i I d D xC 



By using the identities 



T>(ifnf)) e 



d D x{a-g^f 
T>[p) e £ J = cost 

— - / d D x{tx — igifj^^fip) 2 
T>{-k) e 2 J = cost 



(2) 



(3) 



and by inserting them in eq.(Q), we can write 



V{ipip)V{o-)V{T?) e 



i I d D x 



1 



if) (id — M + go + ig^^r ■ Tr)tf) — -{a +7?) 



(4) 



By carrying out the integral over the fermion fields, we obtain an effective action in 
which only the bosonic degrees of freedom do appear 



V(o)Vm e iSB = e iW 



By redefining the o field with the following shift 

go — » go + M 



(5) 



(6) 



we can write 

S B = I d D x 



(a 2 + n 2 )(l + 5Z D )- — 

g 



i\ogdet(id + go + ig^r ■ n) (7) 



where 8Zd is an counterterm coming from the renormalization of the fermion loop [0]. In 
NJL 5Z 4 = 0, since the theory is non-renormalizable. In this case the regularization of 
infinities is obtained by means of a momentum cut-off. 
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To simplify the notations we define 



N = N f N c ; X = Ng 2 ; a = 



M 
W 2 



(9) 



By using the formal identity log det = tr log, and carrying out explicitly the trace over 
colour indices, we can write 

N 



S 



B 



d°X 



2A 



\l + 5Z D ) -NafiSjx 



— iN c Tr \og(id + (f) ■ a) 



= J d D xC B ((f)) 



(10) 



At this point, by following the method of ref.||14j|, we couple an external constant source 
J to the field <p(x), and expand <p(x) around its vacuum expectation value (p 



x 



[x) + ip 



$[J] = (O + |0(x)|O-). 7 



6W 

u 



(11) 

(12) 



Since J is constant, from the Lorentz invariance of the vacuum, it follows that (p is 
constant too. By supposing <f[J] to be invertible, we can define the effective action T[0\ 
as the Legendre transform of W [J] 



n 



W[J{0\\ - I d D xJ[(p](p = -V{(p) / d u x 



(13) 



It can be easily shown that, at vanishing sources, <p must satisfy the stationary equation 



8(p 







(14) 



By using eq.(p^), W[J] can be written as 

W[J] = W [J]+Wi[J] 

where 

W [J] = j d D x(c B (^) + J^) 
and Wi[J] satisfies the integral equation |I4| 



Wi[J] = -zlog J P(0')exp i J, 



i I d D x I C 



<p)-<jt 



, SWx 
6$ 



(15) 
(16) 

(17) 



and C§ is obtained from the lagrangian Cb{4>' + V 9 ) after having subtracted constant and 
linear terms in 6' . 
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By using eq. ([□]), eq.(]TO|) becomes 

N 



S B 



We now define 



d D x 



- — (ft + ^) 2 (1 + 5Z D ) - Na(ft + ipy&ji 



-ziV c Tr \og(id + (p • a + ft ■ 



(19) 



(id + (p ■ a)G(xi — x 2 ) = 5 D (xi — x 2 ) 
The last term in eq.fll8|) becomes 

- iiV c Tr \og(id + tp ■ a + ft ■ a) = -iN c Tr \og(id + <p ■ a) - iN c Tr log (1 + G<f> ■ a) (20) 
Thus 



Wo[J] 



d D x 



^p 2 (l + 5Z D ) - Na&Sjt + Jip 



iN c Trlog(id + (p-a) (21) 



and its final form, at J = 0, is the standard mean field fermionic term 



Wo 



d D x 



N 



~2A^ (1 + 5Zd) ~ Na( f p5 ^ + iN J 



D r d D p 



(27T) 



D 



\og(p 2 + Lp 1 



(22) 



which is the only part which survives in the N c — > +oo limit, since W\ is of order 1/N C 
with respect to Wq. 

To evaluate Wi, one has to expand Sb in eq.(|T8[) in powers of ft, starting from the 
quadratic terms. Now, by ordering the terms of the argument of the exponential in eq. (|T7| ) 
in powers of 1/N C (instead of the h formal series of ref . |[L4| ) , it turns out that the leading 
term is just the quadratic one (this can be verified by the rescaling ft — > ft/y/Nc). 

Thus, by keeping only the leading term of the 1/N C expansion, W\ is given by 



Mi 



= I V(ft) exp {-^ J d D xft 2 (l + 5Z D ) 

f d D Xid D x 2 Tr G{x\ — x 2 )ft(x\) ■ aG(x 2 — Xi)ft(x 2 ) ■ a \ (23) 



By Fourier transforming and evaluating the traces in the previous equation we finally 
obtain 

d D p 



iN 
2A 



(24) 



where j labels each scalar /pseudoscalar field and Dj 1 (p) is given by 
Dj\p) = i[(l + 6Z D )+m j (j>)] 

q-(q + p)±ft 2 



;i + 8Z D ) - iD\ J 



(2tt) d [(p + q) 2 - <p 2 } (q 2 -^ 2 ) 



(25) 
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where ± refer to the scalar and pseudoscalar self energy Hj(p) respectively . 

Notice that, apart from the linear term in (pi, both Wo and W\ depend on the chirally 
invariant combination (p 1 = g 2 (a 2 + vr ), which has to be determined by the eq.Qnp. From 
this condition it turns out that the physical value is for zero pseudoscalar components, i.e. 

g>=(jgv,d) (26) 



3 Free energy density at order 1/N C 

Taking into account the bosonic fluctuations is particularly interesting also as far as 
the thermodynamics is concerned. In fact, although depressed by a 1/N C factor, pion exci- 
tations should dominate the low-temperature equilibrium properties of strong interacting 
matter (this is particularly expected in the chiral limit, although there are limitations in 
applying this method for zero current quark masses in low dimensions, due to well known 
infrared problems). 

The finite temperature extension of the formalism is straightforward. In fact, for 
imaginary times, one has simply to replace 

r r d D r> i n =+°° r ri d r) 

VW-*/A = HW, ; /_IL^_^/_^ (27) 

and thus, at zero sources, the effective potential, defined by the eq. fll3D , is 



V, = -i^ = (28) 

which is the free energy density T . 

By means of the 1/N C expansion discussed in the previous paragraph, JF can be sepa- 
rated in 

J r = J rF + Y j rf (29) 

3 

with 

F F = ^ 2 (l + ^)+^-^E/7^l0g(^+P e +-n) 

= ^-wi^M i+ ^^ 2 ) (3o) 



^ = 27?/(^ l0g [f^ 1(P °"' Pl 



(31) 
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where the inverse bosonic propagator can be calculated at finite temperature by means of 
standard methods 115 



d d q n F {E q 



D- 1 (p 0nl p)=D 1 (p 0n ,p)-tD\ /."" " /lJ - 
.D\ f , o \ f d d q 



+i 



+ 



4 

n F (E q+p ) 



~M 



(2Tl) d 

n F (E q ] 



E n 



(2 



7T 



E„ 



(E q -p 0n y-E, 



2 

q+p 



(E q + POn y - E* q+P 



E, 



q+p 



+ 



(E q+P - p j 2 - m (E q+P + POn y - ei 



where 



and e 2 M 



Po n = iv n = 2nni/p ; n F {E) 



1 



eP E + l 

Aip 2 for the scalar and e 2 M = for the pseudoscalars. 



P 2 + if 2 



(32) 



(33) 



The term T F is the standard mean field term coming from the one-loop calculations 
and it is purely fermionic, whereas the second term, Tb, is the contribution of bosonic 
fluctuations to the free energy Expression (j3lD and ( ^0|) put in eq . (|29D give the equation 
of state when evaluated at the solution of the gap equation fllil) at finite temperature, 
which is 



Sep 







(34) 



In particular, since the free energy F(<-p) is a convex function of (p, (p must be the 
absolute minimum of T and this gives the evolution of the fermion condensate with tem- 
perature. The procedure adopted as well as the expressions obtained are general, although 
the bosonic part of the free energy may be difficult to evaluate in the form of eq. (|31~1) , 
because of the sum over discrete energies po n . It is more convenient for calculations to 
transform the sums into integrals by a two-steps standard procedure (see for instance ||T5fl ). 



The first step is to use the residue theorem to trade the bosonic sums for an integral over 
a circuit around the imaginary axis, and the second step is rotating along a circuit around 
the real axis. 

Namely one can use the relations (ns(z) = l/[exp(/3z) — 1]) 



P n 



lim 

e^o 2ni 

\ r+oo 

lim — / du ubUo) [f(u + ie) — f(ui — ie)] + 



dzf(z)n B (z) + / dzf(z)n B (z) 



2tt 



dujf{iuj) (35) 



provided f(z) exp(— ft\z\) vanishes sufficiently fast at infinity, f(z) has no singularities on 
the complex plane apart from a cut along the real axis, and f(z) = f(—z) (as it turns out 
to be the case for the inverse propagator iD^i^z^p)). 
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Therefore, the bosonic term can be cast in the form 



T 



B 



lim 



d d p 
(2n)° 



+ - 



d d p 

2tt J (2n) a 



duo n B (uj) log 
duo log 



iDj l (uo + ie,p) 
iDj l (uo — ie,p}_ 



iD- 



iuo 



(36) 



apart from infinities independent on j3 and ip. 
Furthermore, we notice that the function g(z) 



iD l {z,p) satisfies the Schwarz re- 



flection principle, g*(z) 
P, we can write 



g(z*). Thus, by reintroducing the explicit dependence on ip and 



, 1 
lim — 

1 



d d p 

(27l) d 

d d p 
(2vr) d 



duo n B (uo) farg iD^iu + ie,p;(p, (5) — n^j 

■oo r -, 

(iculog iDj l (iuo,p\ip, (3) 



(37) 



where arg[/] is the argument G [0, 27r) of the complex number / = |/| exp(i9). 

We notice that (i) the second term of eq.(|37D would be the zero temperature term for 
a free boson gas, whereas now it is temperature dependent since iD -1 is a function of /3; 
(zz) this term may need renormalization at T = 0, due to momentum divergencies [1TB 
(in d = 3 divergent integrals are understood to be regulated by a momentum cutoff). The 
first term, instead, is finite at any T, and tends to zero for T — » 0. 



4 Conclusions 

We have presented a functional integral procedure to calculate the contribution of 
bosonic fluctuations to the partition function in theories with four-fermion interactions 
by means of a 1/N C expansion. Such a procedure has the advantage, as compared to 
diagrammatic methods, to automatically guarantee conservation of the system symmetries, 
through the consistent use of the effective action formalism. 

Our expression for the bosonic term in the free energy is given in eq.(37), as a sum of 
a finite term which vanishes at zero temperature and and of a term which does not vanish 
in that limit, where it must be regularized. 

Fluctuations around the mean field are of greatest physical importance in connection 
to chiral symmetry breaking, manifesting themselves through bosonic excitations, some 
of which, by virtue of the Goldstone theorem, correspond to the lowest bosonic mass 
excitations, and as such they dominate the thermodynamic behaviour at low temperatures 
where chiral symmetry is still far from being restored. 
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Note added: After completing this work we have seen a recent preprint by E.N.Nikolov 
et al. [|hep-ph/9602274| where the functional formalism is applied to the NJL model at 
zero temperature. 
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